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Abstract. We examine the scaling behavior of the entanglement entropy for the 
2D quantum dimer model (QDM) at criticality and derive the universal finite sub- 
leading correction Jqcp- We compute the value of Jqcp without approximation 
working directly with the wave function of a generalized 2D QDM at the Rokhsar- 
Kivelson QCP in the continuum limit. Using the replica approach, we construct 
the conformal boundary state corresponding to the cyclic identification of n- 
copies along the boundary of the observed region. We find that the universal 
finite term is "fQCP = In il — 1/2 where R is the compactification radius of 
the bose field theory quantum Lifshitz model, the effective field theory of the 
2D QDM at quantum criticality. We also demonstrated that the entanglement 
spectrum of the critical wave function on a large but finite region is described by 
the characters of the underlying conformal field theory. It is shown that this is 
formally related to the problems of quantum Brownian motion on n-dimensional 
lattices or equivalcntly a system of strings interacting with a brane containing a 
background electromagnetic field and can be written as an expectation value of a 
vertex operator. 
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1. Introduction 

Entanglement is one of the strangest features of quantum mechanics. Few would 
disagree that it distinguishes quantum mechanics unequivocally from classical physics. 
While it has played an essential role in the understanding of quantum mechanics, 
paradoxically it has been mostly absent from the theory of quantum phase transitions. 
As in the case of classical critical phenomena, the theory of a quantum phase transition 
is naturally based on the description of the scaling behavior of local observables. 
Indeed, much of the theory of quantum criticality is based on this extension of 
the framework of classical criticality allowing for the natural dynamics specified by 
the quantum mechanical Hamiltonian of the system. [T] In this standard description 
the main effect of quantum mechanics is captured by the dynamic critical exponent 
associated z that governs the relative scaling of space and time and is determined 
by the dynamics of the quantum system. It is a major theoretical challenge is 
to understand the role of observables that have no classical analogue in quantum 
criticality. Quantum entanglement has been proposed as a candidate for such a 
measure. [2] 

While there are other measures of entanglement, the entanglement entropy has 
proven to be the most tractable analytically and simple to define. One begins 
with a pure state of a bipartite quantum system A U B. By restricting oneself to 
observing either A or B subregion, the subsystem is described by a mixed state 
with a non-trivial reduced density matrix. Suppose that A is the observed region 
so that the degrees of freedom in B are traced out. The reduced density matrix 
is pa = Trs paub- Non-local quantum correlations are then encoded in the von 
Neumann entropy, S = — Tr pa In pa- For the case of a total system AU B in a system 
in a pure state, the entanglement entropy is symmetric Sa = Sb so that S should 
only depend on common properties of the two regions. Given the non-local character 
of the entanglement entropy, its behavior in extended systems and quantum fields 
theories has a strong geometric flavor. An important early was the realization that in 
systems away from quantum criticality in D space dimensions, systems dominated by 
short distance physics, the von Neumann entanglement entropy scales as the area of 
the observed region, S ~ pi ^ 1 + . . ., where I is the linear size of the region (say ^4) 
and p is a non- universal constant. [3l H] 

The scaling behavior of the entanglement entropy has been studied extensively 
in ID quantum critical systems and it is by now reasonably well understood. Such 
systems are effectively relativistic (z = 1) 1 + 1 dimensional conformal field theories 
(CFT). It has been shown that in ID the entanglement entropy of a subsystem A 
of linear size £ of an otherwise infinite system (i.e. of linear size L — > oo) obeys a 
logarithmic scaling law,[5l [6l [7J [8] S ~ | ln(~) + . . ., where c is the central charge of 
the CFT, and a is the short distance cutoff. The growing popularity of entanglement 
entropy study is due to the fact that many universal properties of quantum systems like 
the central charge, excitation spectra or boundary entropy can be extracted from the 
entanglement entropy without the need to specify an observable: all that is needed 
is a consistent partition of the system. In addition, the dynamical entanglement 
entropy, entanglement generated in a quantum quench, has also been studied recently 
for a number of one-dimensional quantum critical systems. 1101 111! I12j Even in 
strongly disordered quantum systems, which do not have a local order parameter 
such as random spin chains at infinite disorder fixed points [T5] as well as at the 
Anderson transition and quantum Hall plateau transition. |14) the scaling behavior 
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of the entanglement entropy has been shown to be a faithful measure of quantum 
criticality. 

However, much less is known about the scaling of entanglement in spatial 
dimensions D > 1, and how it relates to the scaling properties of local observable. 
While the leading scaling behavior of the entanglement entropy is the area law, its 
prefactor is not universal, except in D = 1 where the area law becomes the universal 
logarithmic scaling. For D > 1 the only universal contributions may only arise from 
sub-leading terms (relative to the area law). The study of these universal corrections 
requires more detailed and subtle calculations. In general the situation is not as clear. 
An exception is the special case of topological phases in D = 2, where there is a finite, 
i.e. 0(£ ), universal correction to the area law which is given in terms of the topological 
invariants of the effective topological field theory of the phase. [THl EE HZ] Recent work 
has showed that several quantum critical systems in D = 2 do not have logarithmic 
sub-leading corrections and, instead, have a universal finite term. [TS1 [THl [201 123 122| 

In this paper we reexamine the scaling behavior of entanglement in a special 
class of quantum (multi) critical points in D = 2 with dynamic scaling exponent 
z = 2. These conformal quantum critical points have the special property that 
the amplitudes of field configurations {\<j>)} in their wave functions have local scale- 
invariant weights. 23 Simple examples of such systems are 2D quantum dimer 
models[23] and their generalizations. (2"51 [231 121)1 |2"T] The norm of these ground state 
wavefunctions is thus equivalent to a partition function of a suitable two dimensional 
classical statistical model at criticality. Labeling a configuration of the classical 
statistical model by the field <fi, this can be expressed, 



With the explicit form of the wavefunction, one can compute the entanglement entropy 
by constructing the replicated reduced density matrix, Tr p\. For these models it was 
shown that the entanglement entropy has a universal sub-leading correction, [19j 



where jqc p = In R for the case of the quantum dimer models at their quantum critical 
(Rokhsar-Kivelson) point. Here R is the compactification radius (see below) of the 
coarse-grained height model, dual to the QDM, the quantum Lifshitz model. [23l [28l [29] 
The same scaling behavior has been shown to hold, within the e-expansion, in 
relativistic cf> 4 quantum field theory, the prototype of a quantum phase transition. |21] 
In Ref. [18] it was shown that the computation of the "spectral moment" of the 
reduced density matrix S n = trpj for the case of local scale-invariant wave functions 
is the same as the computation of a ratio in which the numerator is given by a 
classical partition function of the form of Eq.([T]) for n copies of the system which 
are required to agree on the boundary of region A, while in the denominator no such 
constraint is imposed. Posed in this way, the computation of S n , and the von Neumann 
entanglement entropy S = lim„_>i jh^{S n — 1), is a problem in a "replicated" classical 
2D critical systems with a "conformal defect" along some curve, [30] i-e a problem 
in boundary (Euclidean) conformal field theory (BCFT) on a somewhat unusual 
manifold. The scaling behavior of S n with sub-system size is then not as surprising 
as it is strongly reminiscent of the finite-size scaling of the free energy in large finite 
critical systems. [3TJ [32l [33l [34] On the other hand, it is known that the structure of 
the CFT strongly determines its possible allowed BCFTs. Thus, the structure of the 
CFT that underlies these scale invariant wave functions must similarly play a key role 




(1) 



Sqcp = ^ + 1QCP + ■ ■ ■ , 



(2) 
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in the behavior of the entropies S n . This natural relation of the entanglement entropy 
with boundary (or surface) critical behavior on a suitably defined manifold appears 
naturally in the path integral formulation. [7J 

In this work we use BCFT methods to reexamine this problem for the case of the 
quantum Lifshitz wave function studied before in Refs.[lH HH] paying close attention 
to the role of the compactification radius R by constructing explicitly the boundary 
states of the associated BCFT. [35; Here we show that in the limit of a large aspect 
ratio, L ^> I ^> a (where a is the short-distance cutoff and L is the linear size of the 
full system) the entanglement entropy has the universal finite term jqcp — In R, in 
agreement with the previous result of Ref . [19] . 

On the other hand, numerical estimates of the sub-leading corrections of the 
scaling of the entanglement entropy in the 2D QDM, obtained using the wave function 
of the lattice model and using extrapolation methods, suggested that there may 
be additional universal contributions to 7qcp-[M1 HH] These authors attribute the 
apparent disparity to the boundary conditions used in Refs. [181 [T5] . This motivated us 
to reexamine the BCFT of this problem and to give an explicit form of the boundary 
states. In BCFT, boundary conditions are in one to one correspondence with the 
primaries of the bulk CFT.[35] For the case of multiple copies at hand here the bulk 
CFT possesses additional symmetries. The imposition of boundary conditions along 
a specified curve, a "conformal defect," reduces these symmetries. A similar situation 
occurs when considering defect lines in the Ising model. [30] 

A better understanding of systems of n-coupled critical systems is also crucial 
for reasons outside the entanglement entropy. From a broader viewpoint, a familiar 
condensed matter context are defects in lattice models. As mentioned earlier, one 
such example are defect lines in the Ising model. [30] More generally, conformal defects 
are very hard to classify. Even for two copies of a free bosonic theory, central charge 
c = 2, the complete classification remains elusive. [37] Theoretically, the classification 
of such defects is an important issue. These represent fixed points of a BCFT. A 
familiar example is the £;-channel Kondo problem. [38] More complicated examples of 
n intersecting theories also appear in a more applied context. The intersection of 
multiple quantum wires at the same point is one example. [39 As mentioned earlier, 
systems of n-coupled degrees of freedom also make their appearance in the theory of 
quantum Brownian mot ion [37] 00] and in the dissipative Hofstadter model. [4T ] [42 ] [43] 
In this work, we add the entanglement entropy in conformal quantum critical points 
to this list. 

We compute the quantity Tr p n A in terms of the original degrees of freedom and 
show that there is no additional factor as suggested by Ref. [20]. While our result 
agrees with our previous work, we note that this happens only in the asymptotic 
limit of a long cylinder and that for finite sized systems, there is generically a non- 
trivial n-dependence. To arrive at these conclusions, we identify a new boundary state 
describing the boundary condition that n copies of a system are stitched together at an 
interface, and we argue that additional sub-leading corrections to Eq.© are a result 
of differing boundary conditions used. We show that the correct boundary condition 
for the original degrees of freedom is one where the extra factor vanishes. This turns 
out to correspond to a subclass of conformal defects, at the common boundary. In 
addition, we show that there is a geometrical interpretation to Tr p A as the ratio 
of classical partition functions defined on different tori, and that the universal sub- 
leading term is the "thin torus limit" of this ratio. With an explicit calculation, we find 
analytically that in the n — > oo the entanglement spectrum is given by the characters 
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of the underlying conformal field theory describing the ground state wavefunction. 
This provides a case where it is possible to verify analytically the conjecture of Li and 
Haldane.[33] Finally, we demonstrate that in a string theory language, it is possible 
to think of the common boundary condition as a brane with a background gauge field 
and that Tr p A is the expectation value of an appropriately defined vertex operator 
(in the "target space" not in the CFT of the world-sheet). 

The paper is organized as follows: in the next section, we review the construction 
of Tr p\ for conformal quantum critical models. We then specialize to the quantum 
dimer model in Section |2~T1 with the bulk of the calculation in Sec 12.21 A brief review 
of boundary conformal field theory in | Appendix A is given as it plays a central role in 



the main calculation. Using the explicit expression for Tr p n A we show in Section 13.1 
that the entanglement spectrum has a level degeneracy indicative of the underlying 
conformal field theory describing the ground state wavefunction. In Section 13.21 we 
relate the boundary condition to a brane with a specific magnetic and electric fields 
and show that Tr p A can be computed as the expectation value of a vertex operator. 

Note added to the text: While this paper was being refereed, we received a 
preprint by Oshikawa [45^ and became aware of a number of inconsistencies in our 
treatment of the boundary state. In particular, it was pointed out that our glueing 
condition violated a permutation symmetry in the problem. These corrections and its 
consequences are noted in Section [4] 

2. Conformal Quantum Critical Points 

We are interested in the von Neumann entropy for systems at conformal quantum 
critical points. These are systems where the norm of the ground state wavefunction 
is equivalent to a partition function of a two dimensional classical statistical model at 
criticality, i.e. Eqn. ([1]). With the explicit form of the wavefunction, one can compute 
the entanglement entropy. This is the von Neumann entropy of the reduced density 
matrix 

d 

S en t = -Tr p A \np A = ~—Trp A . (3) 

Since the ground state wave function is a local function of the field </>(x), a general 
matrix element of the reduced density matrix is a trace of the density matrix of the 
pure state ^gsM over the degrees of freedom of the "unobserved" region B, denoted 
by 4> B (x). Hence the matrix elements of pa take the form 

^t\PA\4>t+x) = ^J [Dft 

(4) 

where the degrees of freedom satisfy the boundary condition at the common boundary 

r : 

0f|r = ^|r = ^ i+ ilr- (5) 
This problem can be thought of as a problem in boundary conformal field theory by 
letting A parameterize the boundary interaction between copy i and at the strong 
coupling limit A — > oo the copies are required to have the same configuration on the 
boundary and at A — > the copies do not interact at the boundary. Formally, the 
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trace over n copies of the reduced density matrix can be written as the ratio of the 
partition functions in these two limits, 

Trp'X = f^ff (6) 

The theoretical challenge is to compute this ratio. Various formal mathematical 
devices have been devised, but until now a direct approach has been lacking. It 
is desirable to understand the boundary condition in terms of the original degrees 
of freedom since those are most directly related to a physical dimer covering of the 
lattice. 



2.1. Quantum Dimer Model 

Here we consider the simple case where the action in ([T]) is the Gaussian free field 
theory, described by a bosonic field with the property that it is identified on a circle 
of radius R, (p ~ ip + 2irR, 

S=^- ( d 2 x dtpdtp. (7) 



One can think of this ground state (Q} as a superposition lattice configurations in a 
statistical model such as the dimers on a square lattice [24] and its generalizations 
[25] [46] or as coverings in the eight- vertex model [23] with special choice of the Baxter 
weight [47] . 

One can think of this action as describing a 1 + 1 dimensional system defined 
on some manifold, a cylinder here. To connect with our pervious work, we let the 
circumference be of length i and the length L. Periodic boundary conditions are 
imposed in the ^-direction. The field <p(x, t) is a linear combination of the holomorphic 
4>(x,t) and anti-holomorphic 4>{x, t) parts, ip(x,t) — cf>(x,t)-\-(f>(x,t). Holomorphic and 
anti-holormophic mode expansions can be written as 



<j){x,t) = ip + Y n o(t + ix) + ^a k e 2 * k (t+t 

k^O 



- . . _ 2ir _ sr-^ i j 

<p(x, t) = ip + —^oit - ix) + -a k e 

fc^o 



(8) 



where k € 1 and the zero modes are given by 

/ m wR\ ( m wR\ 

^ = {r + ^)^={r-^)- (9) 

The constants tpa are canonically conjugate to these zero modes, [7To,¥>o] = i- The 
primaries of the boson field theory are labelled by the value of the zero modes, (ttq, ttq). 

At the ends of the cylinder, Dirichlet boundary conditions can be chosen tp(x, t = 
±L/2) — 0. On the oscillator modes, one finds that 

a k = ~a{q k . (10) 

Here, q = e 27rlT where r = i^. For the zero modes, Dirichlet boundary conditions set 
the winding modes w — 0. At the boundary t = continuity of the fields implies that 
lim a ^o </?( a ) = <p{— a )- This restricts the modes to obey 

a t = a k ; vo=^o ; ^+ = 770, (ii) 
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X 

Figure 1. We are interested in the long cylinder limit where L > I where L is 
the length of the cylinder and I is the length of the circumference. 

which amounts to simple continuity of the field at the boundary in agreement with 
iJSJl. Now, we choose to observe the lower half of the cylinder (t < 0), region B, and 
compute the reduced density matrix for the remaining region A by the replica trick. 
The fields at the boundary should be continuous with each other as described in ([5]) 
and follow an explicit relationship among the modes given by Computing the 

reduced density matrix is hence a computation of a ratio of partition functions ([6]). 
This is most simply done in the boundary conformal field theory framework. This 
formalism is briefly reviewed in |Appendix A| As it will play a central role, those 
unacquainted should proceed there. More extensive reviews of the subject matter can 
be found in other canonical texts. [15] 

2.2. Boundary States 

To make use of boundary conformal field theory, we fold the system at the boundary 
at t = (see Figure [1]) so that there are 2n cylinders of half the total length L. 
At t = L/2, Dirichlct boundary conditions were imposed on the ends which relates 
the holomorphic and anti-holomorphic modes. The boundary condition (|10p can be 
regarded as an eigenvalue expression for the annihilation operators, and the Ishibashi 
states can be written as the coherent state 

\m)) = J[ exp (ajjQyof) |m; 0), (12) 
fc=i 

where and a l k \m;0) — and the state |m;0) labels states in the Fock space where 
the winding mode w — 0, as required by Dirichlet boundary conditions (fpQ|). With 
two copies, the matrix is given by, Q,j = —q k 5ij. The boundary state is a linear 
combination of these Ishibashi states, given by (|A.6|) and (|A.11I) 

oo 

\B D )= 9D J2 e^°|m». (13) 

m— — oo 

<7d is the ^-factor associated with the Dirichlet boundary condition for the free boson 
and can be computed explicitly as mentioned in (|A.12[) . 

We now construct the non-trivial boundary state that is at the boundary t = 0. 
Below, we specialize to the case of n — 2 copies, but the result is easily generalized 



Universal Behavior of Entanglement in 2D Quantum Critical Dimer Models 8 

to arbitrary n. At A —> 0, the two copies are decoupled, and the boundary condition 
becomes 



it 




2+ 




k 


= a? 




= "fe 






4+ 




r 


= at 




= *2 



(14) 

and the momentum modes are restricted to obey, 













= A 














4 


= 4 


-o 3 


-4 



(15) 

The equation (fT4"|) can be regarded as eigenvalue equations for the destruction 
operators and the Ishibashi state for each can be written succinctly as 

x((m'-w'\ = (m';w'\ JJ exp (4^y(A)a J fc ) , (16) 



k=l 



where 



^0100^ 
10 
1 
\ 1 J 



(17) 



A similar construction can be done for the A — >• oo case. In a previous version of 
the paper, IZij (A — > oo) was given incorrectly and was pointed out by Oshikawa [J5] 
to violate a permutation symmetry of the problem. Instead all possible boundary 
conditions must be summed (see Section [4]). The resulting matrix is given by equation 
(HEJ). 

The state (m'; w'\ denote states with eigenvalues obeying the matrix relationship 

TTi = TZijTTj (18) 

where the form of 7£y depends on the value of A. It is given by p7|) for A — > and 
(|48l) for A — > oo. Hence, the boundary state can be written as a linear combination of 
Ishibashi states again, 

(B X \ = x ((m';w'\gxJ2 e ^' ' P0 - ( 19 ) 
where the sum is over modes on the correct compactification lattice. This subtlety 
was neglected in a previous version of the paper, but was pointed out by Oshikawa 
[4"5] . Here, the g- factor will play an important role. We call the g- factor associated 
with the A — > oo limit guv and the g-factor associated with the A — » limit gm. We 
leave the specific evaluation of these normalization factors till later. 

The partition function is then easily evaluated as the overlap of the two boundary 
states. 

Z x (n) = (B x \q 6 \B D ) =g x g Dq - 2 ^ 12 (i P ' \ q A ^\i Po )Z osc , 

(20) 

where the quantity (tp' Q \q H \tp ) is the piece only involving the zero modes and 
Z osc involves the oscillator modes. The respective matrices describing the glueing 
conditions inserted, one readily finds, 

Z osc (X -+ oo) = J] (1 - (qT^ (! + 1 k )~\ 

k>0 

Z osc (\^0) = H(l-(q 2 ) k y n . (21) 

fe>0 
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Figure 2. In the two limits A — > oo (a),(c) and A — > (b),(d), the 
holomorphic (ingoing arrows) and anti-holomorphic (outgoing arrows) are stitched 
in topologically distinct ways. The legs are labelled tpi,tp2, etc. in a 
counterclockwise way. In (a), there is effectively only one bosonic degree of 
freedom, while in (b) they are stitched for form n-independent closed loops 
associated with an independent bosonic degree of freedom. In figure (c),(d) a 
similar diagram can be drawn to compute the g-factors. In (c), A — > oo, there is 
only one non-chiral bosonic degree of freedom while in (d), A — > 0, there are two 
non-chiral bosonic degrees of freedom. Because of the permutation symmetry, one 
must sum over all possible closed, chiral loop configurations. This corresponds 
to a sum over glueing conditions in Section [4] and the glueing matrix given by 
Oshikawa |45] , 



The zero modes require some care. In a previous version of this paper, the winding 
and momentum modes were assumed to be independent, but the glueing matrix forces 
a enforces a condition. Instead one must sum over a lattice 5 as noted by Oshikawa 
[4"5] . This is dual lattice which satisfies the glueing condition on the zero modes, 

7Tj = IZijTTj (22) 

where TZij is given by ([35]). Hence for A — > oo, one finds, 

(^k^°°Vo) = E ( 23 ) 
By a similar set of arguments, for A — > one finds, 




(<£ok H(A IVo) = I >J (q 2 )^) ■ (24) 

\m— — oo / 

The results are easily generalized to arbitrary n. Writing in terms of ^-functions (see 



Appendix B I , the UV and IR limit partition functions are given by, 

Z\^oo{n) = gvvgo 
Z\^a{n) = gmgD 

As a non-trivial check, for Z\^$ we note that 2r is the modular parameter for a 
cylinder of total length I and observe that the term in paranthesis is the C/(l) character; 
it is the partition function for n-decoupled bosons on a cylinder of length £ as it should 
be. In light of this observation, gjn should be equal to 1. 
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The factors guv and gin are normalization factors that have yet to be fixed. 
These can be fixed by Cardy's conditions and one finds that, 

9IR = 1, guv = R {1 - n) Vn~ (27) 

Crucially, the ratio guv/giR = i?^ 1_n ' \fn 1 and hence we arrive at the main result 
of this paper, 

Tr p n A = C(n,£,L)^(^) (28) 

9IR \ ^A^O / 

where Za-j-oo and Z\^o are given by (|25[) and ([26]) respectively and we have factored 
out the constant factors guv/giR explicitly. Here, we have included the regulator 
dependent contribution in the function C(n,£,L) that comes from the short distance 
cutoff which we have hitherto suppressed. 49 We comment on its specific form below. 
More importantly, we note that the each of terms $3/77 are characters of a U(l), c = 1 
conformal field theory. [48] 

Thus far, we have been concerned only with the regulated part of Tr p n A and have 
neglected the divergent contributions from the short distance cutoff. These can be 
recovered by a careful £- function regularization. [15] In general, it was shown that in 
the limit of interest the free energy for a system on a smooth open manifold scales 
as,[33 j 

In Z x = n a \A\ + n h P - C - X In* + 0(1) + 0(£/L). (29) 

Here A is the area of the cylinder, P the perimeter of the boundary and x is the Euler 
character of the manifold (zero for cylinders). The coefficients p, a and depend on 
the short distance cutoff. The order one contributions are what have been computed 
explicitly thus far. We note in both limits, the total area is simply nA where A is 
the area of a single cylinder so that in the ratio ([6]) this contribution cancels. The 
difference in the two limits lies in the perimeter of the boundary. In the decoupled 
limit, there are 2n unshared boundaries of length £ and n shared (where the two 
halves are joined) boundaries so that the total perimeter of the boundary is (n + 2n)£. 
Meanwhile, in the strongly coupled regime all the manifolds coincide (smoothly) on a 
single interface of length £ and there are still the same 2n unshared boundaries. The 
total length of the perimeter in this case is (1 + 2n)£. The important point is that 
there is also a divergent non- universal contribution to Tr p\. To leading order, 

C(n,l,L) = e ^(i-n)e_ (30) 

where the prefactor pi, is non-universal and depends on the short distance 
regulator. [49] It is important to note that this analysis neglects the effects cusps and 
corners that occur if the boundary is not smooth, giving rise to conical singularities, 
but these have been shown to give £ dependent contributions that scale as a power 
law that may be non-trivial functions of n and to additional universal (logarithmic) 
corrections to the entanglement entropy. [HJ [50] 

It is important to observe that as n — > 1, Tr pa = 1 as it should be. Now, in the 
present long cylinder limit the ^-functions and ^-function are equal to the identity 
to leading order in e - 27Tl / T . Our main result is that the universal sub-leading term to 
the entanglement entropy is, 

7QCP = ln-R- -, (31) 

where, once again, R is the compactication radius of the boson. This agrees with 
numerical result for the quantum dimer model. [20] 
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3. Discussion 



While the result for jqcp is the same as our previous result, [19] the interpretation is 
much different. The full analytic expression for Tr p A , for instance, differs and it is 
only in the asymptotic limit that the two results are identical. In the very long cylinder 
limit, L/l — > oo the characters asymptotically approach the identity, and one finds that 
Tr p n A identically reduces to the one obtained previously. [19] However, for large but 
finite values of the aspect ratio L/£, Tr p n A has a non-trivial n-dependence. Naturally, 
an immediate question is what information is contained in Tr p n A in this picture? We 
show below that the entanglement spectrum is indeed described by the underlying 
conformal field theory describing the ground state wavefunction. More importantly, 
several immediate questions arise. First, the form of Tr p A is relatively simple and one 
wonders if there is a deeper reason for this. Secondly, we have constructed boundary 
states |£?a-s.o) and |Z?a->oo) which lack a straightforward classification as free or fixed. 
They represent instead a coupling between different copies. An understanding of these 
states is clearly desirable. Lastly, our result for 7qcp hinged on guv/giR = R^ n ^ \fn 
and a natural question is whether other values are possible. In what follows, we address 
each of these issues. 



3.1. The Entanglement Spectrum 

With an explicit expression for Tr p A , one in fact can construct all the moments of 
the Renyi entropy, 

S n = -^—\nTY P n A . (32) 
1 — n 

One can examine the degeneracy of states by examining the higher moments, n — > oo, 
of the Renyi entropy. It has been postulated, but not shown analytically, that 
the higher moments of the Renyi entropy should be given by the characters of 
the underlying conformal field theory describing the ground state wavefunction. [33] 
The characters of a conformal field theory count the number of independent states 
occurring at a given energy level. Here, we find that in the n — > oo limit the 
contribution from the strong A — > oo coupled sector tends to the identity so that 
the Renyi entropy in this limit is given by, 

= ^— log Tr p\ = -J- (~^ + log Hq)) +... (33) 
1 — n 1 — n \ 12 l J 

where the . . . indicate subleading constant contributions. Here, <j){q) is the Euler 
function that is related to the 77- function through Ramanujan's identity. 

m - r c/3 V?)- (34) 

l/4>(q) is related to the partitions of integers, p(k), and q — er 2lxl l T where r\{— 1/2t) = 

1 OO 

M = X>)? (35) 
In the long cylinder limit, q is exponentially small and 63 —> 1. One finds that, 

S n ^ 00 (L>t) = ^-^j-]n(jr P (k)q k ) . (36) 



\k=0 
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Asymptotically the multiplicities are given by the the partition of integers which 
exactly describes the number of states at a given energy in a c = 1 free bosonic 
conformal held theory. Therefore, in a hnite size computation of the entanglement 
spectrum the degeneracies of the eigenvalues of the reduced density matrix are 
(asymptotically) given by the integers p(k). This relation was conjectured (for 
quantum Hall wave functions) by Li and Haldane. [44] 



3.2. Background Electromagnetic Fields 

One issue still untouched by the preceding discussion is where additional universal 
corrections to ^qcp might come from. Some insight can be gained by realizing that 
the problem of 2n free bosonic held theories interacting only at a common boundary 
has been studied in the context of quantum Brownian motion in a magnetic field. [3 7) 
and open strings in a background gauge field. [HJ [51] The connection can be seen 
more concretely by considering Tr p\ in the path integral formulation 

Tr Pa = ^F-FT JTTT . (37) 

where = li m A->oo an d So[<p] = limx->o Sx[<p] describes the bosonic action 

of n scalar fields satisfying boundary conditions specified by A. Once again, the 
numerator describes n fields that are forced to coincide at the boundary (hence A — > oo) 
and in the denominator the n fields are decoupled from each other (hence A — > 0). 

We will now write S[ip] in a form that we find more useful as follows. The first step 
in understanding this problem in terms of quantum Brownian motion in a magnetic 
field or, equivalcntly, open strings in a background gauge field, is to fold the system 
accross the boundary, thus doubling the number of fields. Let with i — 1, . . . , 2n 
denote a 2n component scalar field whose upper n components label the (folded) fields 
from the A regions and its remaining (lower) n components are those of the B region, 

$ l = (^,...,^,<^,...,^). (38) 

The action for the $ field is (for so far unspecified boundary conditions) 

sm = - d* x j2( d ^f- ( 39 ) 

** i=l 

Here we choose our coordinates so that x is the direction parallel to the length of the 
cylinders and t runs along the circumference. 

Now we make a "T-dual" transformation on the B field, ipg \ . . .tps™ ■ This 
corresponds to a symmetry of the action (on the B fields) with respect to the 
interchange of their winding and charge modes. [48l [51] 

i^ w i> ^ m -< (40) 

Here m 1 and w l label the winding and charge numbers of the zero modes of the 
fields tp l B . Under T-duality the compactification radius R transforms as R — > 2/R, 
while Neumann boundary conditions transform into Dirichlet boundary conditions, 
and viceversa. 

The fields only interact with each other at the common boundary via the boundary 
conditions. To this end we introduce a field Aj = ^Fij& localized at the boundary, 

Sm = lf<f x YK 9 ^) 2 - <fdtAd t <P l . (41) 

i—1 
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where F^ is an antisymmetric matrix we define below. Formally, this action describes 
the "dissipative Hofstadter model" [13] (with vanishing potential). Upon varying the 
action, the fields are found to obey the usual wave equation with the boundary 
condition, 

d x & F,,ih<\>-. (42) 
If Fij is a constant matrix, independent of and anti-symmetric, then clearly 



F 



One can think of F^ as the 2-form field strength tensor associated 



with a gauge field Ai. Now, letting F^ be the 2n x 2n matrix, 



F — 

r% 3 — 







-M 1 




(43) 



it is readily seen that Eq. (|4"2l yields the desired boundary conditions at the common 
boundary for the scalar fields by a suitable choice of the n x n matrix M. 

This construction can be used to represent both the numerator and denominator 
of Eq. ([57)1 by the choices (in this basis) 

/ 1 2 2 
2 1 2 



2 2 1 



\2 ■ ■ 

for A — !• oo and A — > respectively. 



\ 



Mn 



/ 1 





V o 



o \ 





(44) 



Previously, a matrix was given for M^ that 
violated the permutation symmetry present in the boundary condition. Summing over 
possible glueing conditions yields the matrix above. The matrix M m is then identical 
to the matrix 7Z (see equation ([48)) ) of Oshikawa [45] up to an overall unimportant 
sign. 

In analogy with electromagnetism, there is a electric and magnetic field at the 
boundary x — that mediates the interaction between the 2n-dimensional systems. 
In terms of quantum Brownian motion, [37] one can regard the problem as a particle at 
x = on a 2n-dimensional space moving in a electric and magnetic field. As an open 
string, one sees that the boundary can be thought of as a brane carrying a magnetic 
and electric field. |54) [42] The analogy also gives a possible interpretation of the result 
obtained by Ref. [20) . In the context of quantum Brownian motion in n-dimensions, 
it has been shown that there is a plethora of possible boundary states where guv I 9ir 
is non-trivial. These boundary states correspond to a different electromagnetic field 
and hence a different coupling of the replicas in our picture. 

We have thus constructed a gauge field (and associated 2-form) that describe the 
boundary interaction. One can go further by rewriting Tr p n A in the suggestive form 

^PA- \ efd ^>y (45) 

where A°° corresponds to the choice M m and A® to Mo (both defined above). 
The universal sub-leading correction to the entanglement entropy is the asymptotic 
behavior of this correlation function. In string theory, such objects are generically 
called vertex operators. [ST] This operator can be understood as counting a topological 
charge 



(46) 
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where the field is a map, : T 2 — > T 2n . Integrating by parts yields the correct 
boundary field. Importantly, there exists a basis of the field where is an n- 
block diagonal matrix of anti-symmetric tensors. $ J can then be written as a 
tensor product of maps T 2 —> T 2 and Q can be thought of as a product of holonomies 
characterized by the homotopy group 71^2 (T 2 ) = Z. An interesting observation is 
that the gauge field Ai is fixed by the free part of the action of Eq. (|39|) so that a 
perturbation that brings the system into a topological phase, e.g. mcos(ip), has no 
effect on the boundary condition and hence Ai remains the same. We expect that this 
correlation function, and hence the entanglement entropy, should take different values 
in the different phases. 



4. Note Added to the Text 



While this paper was being put into production, we received a preprint by Oshikawa 
|45] and we became aware of a number of inconsistencies in our previous work. In 
particular, the matrix used to described the coupled boundary condition did not 
possess the full permutation symmetry of the problem and the winding and momentum 
quantum numbers of the zero modes were not treated correctly. As noted by Oshikawa, 
|45j there is a permutation symmetry in the problem and each copy and definition 
of region A and B are interchangeable. One should then sum over all possible 
combinations of matrices Pij(X — > oo), where Pij(X — > oo) is a shift matrix. The 
glueing condition is then, 

Kij = l + (P l] +P} j ) + P 2 ] +P^ + .. 

The resulting matrix 1Z is in fact the matrix 1Z given by Oshikawa 
unimportant sign. Explicitly, 



p(2n-X) p (2 

ij ij 



-1) t 



(47) 

up to an overall 



Hi 



( 1 


2 2 . 


• 2 \ 


2 


1 2 . 


. 2 


2 


2 1 . 


. 2 






. 2 


V 2 




• 1 / 



(48) 



Using the sum of all possible glueing conditions because of the permutation 
symmetry, the partition function (|25[) is given by the expression found by Oshikawa 
145]. 



Z\^>oo — 9uv9d 



1 \ i /o/i "T~r 1 



--1/24 



TT ^—Y 

A J- 1 + q m 



q K2 / 2 (49) 



m=1 " Kes 

Here, 5 is a lattice that satisfies the boundary condition on the zero modes, 

7Tj = TZijTTj (50) 

The result for Z\^ is unchanged and given by equation (|2"o]l . As pointed out by 
Oshikawa 1451. the correction normalization is, 



g uv = R 1 - 71 ^ 



(51) 



The universal subleading term to the entanglement entropy is dependent on this ratio 
of g-factors so that the result (f5T|) is now 

S = logR-l (52) 
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Several results in the paper are modified because of this. In particular, the 
quantity Tr p n A (i.e. equation (|28p ) is now given by 

_J_\ ~— 1/24 TT 00 1 %K 2 /2 

X1 v a — ; r^n l 5 "^ 



Then in the limit L 3> £, q — > one finds that the distribution of eigenvalues of the 
Renyi entropy for large moments n, is given by 

S^oo - ^ ^ Tr ft = T ^(-~ + ^ m) + (54) 

where 0(g) is the partition of integers and . . . include subleading constant terms. 
Eigenvalues are distributed according the level counting for the free c = 1 boson 
theory that describes the ground state wavefunction, as mentioned in the body of the 
paper. 



5. Conclusion 



Before concluding with broader and more speculative issues, we summarize our results 
for the quantum dimer model. We found that in the limit L 3> £, Tr p A is given by the 
expression of Eq. (|53j) . This was done in terms of the original degrees of freedom and 
it is hoped that this clarifies what boundary condition must be used at the common 
interface. In the limit n — > 1 of S = — <9„Tr p\, it was found that the universal finite 
part of the entanglement entropy is 

7QCP = lni?-i. (55) 

This result coincides with the numerical result. [5D] We found that for finite sized 
systems, there is a non-trivial n-dcpcndcncc, not reflected in the construction of 
Fradkin and Moore. [TS] The source of this difference laid in the subtle details of 
defect lines in critical systems. In this work, we showed that the boundary condition 
is not described by the same notion of Dirichlet boundary conditions as in the original 
system. We circumvented the difficulties here by working directly with the original 
degrees of freedom. 

By considering the original degrees of freedom, we further demonstrated that in 
the limit n — > oo , Tr p\ has a distribution of eigenvalues characterized by the correct 
underlying conformal field theory of the ground state wavefunction, confirming, at 
least for this case, a conjecture put forth by Li and Haldane.pH] Attempting to 
understand where the universal corrections to entanglement entropy come from, we 
related the problem formally to work done on quantum Brownian motion and branes 
with a background electromagnetic field, and showed that Tr p n A can be understood 
as an expectation value of a vertex operator. 

In this work we focused on the quantum dimer model, i.e. models where the norm 
of the ground state wavefunction is related to the free Gaussian field theory, but the 
methods can be readily extended to different conformal quantum critical models. It 
would be interesting to see if a similar structure exists for more complicated systems, 
such as non trivial topological theories with non-Abelian excitations. [5S1 [SHI HH1 HZ] 
For some simple cases, an exact solution is possible. [57] 
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Appendix A. Boundary Conformal Field Theory 

To exploit conformal invariance, it is useful to think of the system as being on a 
cylinder with circumference j3 and length I with boundary conditions A, B on the 
field on the left and right ends of the system respectively. Quantum mechanically, this 
corresponds to evolving a one dimensional system defined on the line x in time /3. The 
partition function is given by the usual expression, 

Z AB = Tie~^. (A.l) 

If the Hamiltonian possess conformal invariance, then one knows that time and 
space can be interchanged, iflior equivalently, the system is invariant under the 
modular transformation S. One now has a cylinder which is wrapped around in the 
spatial direction and extending upward in time. The corresponding Hamiltonian in 
this picture can be regarded as propagating the system for the time interval I from 
the initial and final state \A), \ B) 

Z AB = (A\e- eH "\B). (A.2) 

The states here belong to the Hilbert space of states quantized on the circle, i.e. they 
can be decomposed into linear combinations of states in the representation of the 
Virasoro algebra which are labelled by (h, h), the highest weights. [48] 

Because conformal invariance is so restrictive in two dimensions, one can say 
more about the boundary states \A). One typically imposes the condition that 
T(z) — T{z) where T, T are the holomorphic and anti-holomorphic components of 
the stress energy tensor, z = t + ix. In the x, t basis, this means that the diagonal 
components of the stress energy tensor vanish at the boundary T x ^. If the boundary 
is in the time direction, this means no momentum flows out of the system. The stress 
energy generates the conformal symmetry so that the boundary states must satisfy 
the condition, 

[T - f] \A) = 0. (A.3) 

Fourier transforming, this can be written in terms of Virasoro generators, 

[L n - L_ n ] \A) = 0. (A.4) 

This implies that the boundary state \A) must made out of states with the holomorphic 
and anti-holomorphic sectors stitched together in a specific way, i.e. 



\h)) = 2^\h;m)<S\h;-m). (A.5) 



Here, m labels the descendant level in the representation h that belongs to the subset of 
representations that appear simultaneously in the holomorphic and anti-holomorphic 
sectors of the Virasoro algebra. The state \h)) are known as the Ishibashi states. [55] 
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The Ishibashi states turn out to form a basis for the possible boundary states, and 
one can write an arbitrary state \A) as a linear combination of the Ishibashi states, 

\A)=J2cm- (A.6) 

i 

Hence, the characterization of a conformal boundary condition is reduced to finding 
the matrix elements C\. Now, using the expression (IA.6I) into (|A.2[) one finds that 

Z AB = Y d C i A Cf{{i\e- lHf, \€ ) ). (A.7) 

i 

The overlap can be identified with the character of the representation i, \i (e -4 *'^/' 3 ) . 
Now noting that (IA.1|) can be written as a sum of characters, and using the fact that 
the two quantities are in fact equivalent by conformal invariance leads one to the 
relationship 

£ * (e- 4 ^) = £ n ABXl (e-«/') , (A.8) 

i i 

where n AB are the multiplicities that indicate the number of times a representation i 
appears in the Hilbert space with boundary conditions A and B. One then notes that 
the characters transform among themselves by the modular S'-matrix so that 

Xi(e-^ e )=j:si Xj (e-^). (A.9) 

3 

If the characters are linearly independent, then this leads one to Cardy's equation 
which relates the multiplicities n AB that characterize the spectrum of the theory for 
fixed boundary conditions A, B and the matrix elements C\ that characterize the 
boundary states, 

Y. S i ni AB = C A Cf. (A.10) 

3 

The key problem in boundary conformal field theory is finding a set of boundary states 
where the multiplicities are non- negative integers. 135] [59] For the free boson, a solution 
to this requirement is 

oo 

C A oc e l *°^ V0 , (A.ll) 

w,n— — oo 

where ir and (fo are defined in ©. 

Note that linear combinations of boundary states \A) also satisfy the above 
constraints. An additional choice that is imposed is that n° AA — 1, that is to say that 
the identity representation appears exactly once in the spectrum of the theory with 
A, A boundary conditions. Operatively, this fixes the normalization of the boundary 
states so that in the long cylinder limit Z AA contains the identity exactly once. This 
gives the ^-factor of the boundary state. [60] 

g A = (0\A). (A.12) 

It has been conjectured that relevant boundary perturbations drive the system to fixed 
points given by lower values of the g- factor. [50J In this sense, the g- factor is also a 
characteristic of the boundary condition. 
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Appendix B. ^-functions 
The ?9-functions are defined as 

oo 

0i(i/|t)=i {-l) n q^ n - 1/2)2 e™ (2n - 1)u (B.l) 

n— — OO 

= 2a 1 1 ' 8 sin(7ri/) /(g) Y[ (1 - 2cos(27ri/) g m + q 2m ) , 



m—1 

\2 • 



M"\r)= Y, g 2(«-l/2) 2 e -(2™-l), (B2) 

TL — — OO 

OO 

= 2qV 8 cos{irv) f(q) TJ (l + 2 cos(27r% m + g 2m ) , 

m—1 

oo 

<?3(^|r) - 51 ^"V 21 ™" (B.3) 

n— — co 

OO 

= /(<?) I] ( X + 2 COs(2^ J ,) g "- 1 /2 + g 2n-l^ ^ 
m—1 

oo 

^4(^|r) = J2 {-l) n q? n2 e l2mv (B.4) 

n— — oo 

OO 

= /(<?) II ( : ^ 2 cos(2^ i ,)g"- 1 / 2 + g 2 "- 1 ) , 

m — 1 

where g = e 2lTlT and 

«-np-fl-(^Ti-r- (B - 5) 

m—1 

The r/(r) function is then denned as 

V(t) = q 1/M f(q). (B.6) 

The action of the modular transformation 5:r^-l/ron ^-functions can be found 
by making use of the Possion resummation formula, 

OO _| OO 

E -Trn 2 A+2mrAs _ ■kAs 2 -ttA^ 1 m 2 -2iixms 

TL — — oo v m — — OO 

(B.7) 
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